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Abstract 

Simultaneous multifrequency observations of the Crab pulsar giant pulses (GPs) 
were performed with the 64-m Kalyazin radio telescope at four frequencies u = 0.6, 
1.4, 2.2 and 8.3 GHz, using the K5 VLBI recording terminal. The K5 terminal pro- 
vided continuous recording in 16 4-MHz wide frequency channels distributed over 4 
frequency bands. Several thousand GPs were detected during about 6 hours of obser- 
vations in two successive days in July 2005. Radio spectra of single GPs were analysed 
at separate frequencies and over the whole frequency range. These spectra mani- 
fest notable modulation over frequency ranges, Au, both on large {Au/u ^ 0.5) and 
small {Au/u ^ 0.01) frequency scales. Cross-correlation analysis of GPs at 2.2 GHz 
showed that their pulse shapes can be interpreted as an ensemble of unresolved bursts 
grouped together at time scales of ~ 1 yus being well-correlated over a 60-MHz band. 
The corresponding CP cross-correlation functions do not obey the predictions of the 
amplitude-modulated noise model of Rickett (1975), thus indicating that unresolved 
components represent a small number of elementary emitters. 
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1. Introduction 

Giant pulses (GPs) represent an extraordinary phenomenon in pulsar radio emission. 
First, they have huge peak flux densities, that may be as great as several millions of janskys 
(Cordes et al. 2004; Soglasnov 2006). Second, their duration is ultra short with CP fundamental 
components being even as short as < 2 ns (Hankins et al. 2003; Hankins & Eilek 2007). Third, 
these elementary nanopulses are strikingly almost completely circularly polarized, either left- 
or right-handed (Cognard et al. 1996; Popov et al. 2004; Hankins et al. 2003). Giant pulses are 
reliably known only from seven pulsars (see, e.g.. Knight 2006, and references therein), but the 
properties of GPs were studied extensively only for two objects: for the original millisecond 
pulsar B1937+21 (Cognard et al. 1996; Kinkhabwala & Thorsett 2000; Soglasnov et al. 2004; 
Kondratiev et al. 2006), and for the Crab pulsar (see the most recent review by Bilous at al. 
2008 and references therein). In contrast to regular radio pulses, which typically have Gaussian 
distribution, the cumulative distribution of CP energies follows a power law N{E > Eq) oc 
Eq" with the exponent a being in the range of 1.5-2.5 (Popov & Stappers 2007; Bhat et 
al. 2008). These authors also reported about the deviation from the single power law fit of 
cumulative distributions at low energies. No deviations, however, were found at the highest 
energies even during 160-hours monitoring of the Crab pulsar with the Kalyazin 64-m radio 
telescope at frequencies of 0.6, 1.65 and 4.85 GHz (Popov et al. 2008). These three-frequency 
observations demonstrated that the pulse energy, in general, grows with decreasing frequency, 
and the spectral index is of about —1.8. But radio spectra of individual giant pulses manifest 
deep modulations both at small and large frequency scales (Hankins & Eilek 2007; Popov & 
Stappers 2004; Popov et al. 2008). In this paper we present an analysis of results obtained 
from simultaneous 4-frequency observations of GPs from the Crab pulsar conducted with the 
64-m Kalyazin radio telescope at the frequencies of 0.6, 1.4, 2.2, and 8.3 GHz using the K5 
VLBI recording terminal. The description of the observing setup and data reduction pipeline 
is given in Section 2 and 3. The results of our CP studies at different frequencies are reported 
in Sections 4, 5, and 6, and the conclusions are given in Section 7. 

2. Observations 

Multifrequency observations of the Crab pulsar were conducted on July 20-21, 2005 with 
the 64-m Kalyazin radio telescope simultaneously at four frequency bands: 0.6, 1.4, 2.2, and 
8.3 GHz. The total observing time was about 6 hours. The simultaneity of observations at four 
frequencies was possible owing to the multifrequency feed that the Kalyazin radio telescope is 
equipped with. To record the data we used the K5/VSSP VLBI recording terminal (Kondo 
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et al. 2002) and the K4 baseband converter provided by the Kashima VLBI group of the 
National Institute of Information and Communications Technology (NICT, Japan) for joint 
VLBI observations of pulsars under mutual agreement. The block diagram of the antenna feed 
line with the K5/VSSP VLBI recording terminal is presented in Figure 1. The terminal consists 
of 4 workstations with K5/VSSP PCI-bus boards, that allow to record up to 16 channels, 
each with the bandwidth of 4 MHz. The channels may be distributed in frequencies and/or 
polarizations by using the K4 baseband converter. In our observations we used the frequency 
setup presented in Table 1. the channel frequencies /1-/4 from the Table correspond to the lower 
edge of each frequency channel, and only upper sub-bands were recorded. At the frequency 
band of 600 MHz only one 4-MHz channel was used because of limitations by radio frequency 
interference. In all other bands we used four channels distributed along the whole band as 
indicated in Table 1. 
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Fig. 1. Block diagram of antenna feed line of the 64-m Kalyazin radio telescope with the K5/VSSP VLBI 
recording terminal. 



The system equivalent flux density (SEFD) for every band is dominated by the contri- 
bution from the Crab nebula and is given in the Table 1. Flux densities for the Crab nebula 
were calculated by using the relation Sf = 955/"°'^^, with / expressed in GHz (Cordes et al. 
2004). At 8.3 GHz solid angle of the radio telescope beam covers only 20% of the area occupied 
by the nebula. Thus, we estimated the contribution from the Crab nebula to be 5 times less at 
this frequency assuming an uniform surface brightness distribution. 
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Table 1. Observations summary. Columns list the frequency band, the frequency channels setup for every band with 
the values corresponded to the lower edge of the channel, /1-/4, the system equivalent flux density (SEFD) including the 
contribution from the Crab nebula (see text), the estimated dispersion smearing time in a single frequency channel, ttjmi 
the total number of events found in a total band, A'^p, the peak flux density of the strongest event, Smax- 
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* Giant pulses can be simultaneously detected in more than one frequency channel, 
i.e. the actual number of unique pulses detected is smaller than the total number of 
events in all frequency channels, A^p. 
** See section 6 for explanation. 



3. Data reduction 

The core of the processing pipeline consisted of successive decoding of the data, correc- 
tion for bit-statistics and dedispersion. Every single frequency channel was decoded according 
to the K5/VSSP (IP-VLBI) Board Data format^ to produce the raw binary data stream. The 
decoded data represented the stream of only four alternating numbers: —3, —1, +1, and +3 
according to the VLBA coding of 2-bit digitized values. These data were corrected then for 
the bit-statistics, i.e. the decoded values were adjusted to real ones that correspond to instan- 
taneous ±1(7 brightness levels by using the technique developed by Jenet & Anderson (1998). 
To coherenetly dedisperse the data we used the predetection dispersion removal technique pio- 
neered by Hankins (1971, see also Hankins & Rickett, 1975). More details about the correction 
for bit-statistics and coherent dedispersion can be found in Popov et al. (2002). However, for 
the coherent dedispersion, unlike the Eq. 1 from Popov et al. (2002), we have used the exact 
equations for phase correction: 

,,,,, MM/A/y 

where DM is the dispersion measure in pc cm~'^, D = 2.41 x 10"^^ pc cm~^ s the dispersion 
constant, /o the lowest frequency of the frequency channel, / the current frequency, and A/ = 
/-/o. 

To allow for the direct comparison of times of arrival (TOA) of giant pulses, we reduced 
TOAs in all frequency channels within a band to the highest frequency (except for the data 
at 0.6 GHz where we had only one frequency channel). This was done by introducing a time 
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delay via a linear phase correction kj in the corresponding spectra simultaneously with the 
correction for 5(j){f). The linear correction kj given by kf = 27r x dt/N, where dt is the 
time delay between corresponding frequencies in samples (including fractional part), and is 
equal to the number of samples in the array used for the Fast Fourier Transform (FFT). We 
used = 2097152, corresponding to the time interval of 262.144 ms and covering about 7.8 
periods of the Crab pulsar. The dedispersion correction of every such chunk of data resulted in 
a number of spoiled samples in the restored signal due to the cycle nature of the convolution 
made through the FFT. The duration of this spoiled portion is equal to the dispersion time in a 
single frequency channel, tdm? of which the estimated approximate values are given in Table 1. 
To reconstruct the signal without gaps the successive chuncks of data were overlapped by tdm- 
The used value of DM 56.738cm~^pc was taken from the Jodrell Bank Crab Pulsar Monthly 
Ephemeris^ (Lyne et al. 2008) for the epoch close to our observations. 

To search for CPs at all frequencies except 600 MHz^, we set up the detection threshold 
at a level of 20o" above the mean value for the signal with the original sample tome of 125 ns. 
The statistics of amplitudes x of such a signal is close to a x^-distribution with 2 degrees 
of freedom P{x > a) = exp(— a), where P is the probability of the amplitude x of the signal 
exceeding the threshold a, and a is in units of a. Thus, the probability of false detection with 
a threshold of of 20cr above the mean level is 7.6 x 10~^°, and for the sampling rate of 8 MHz 
it will lead to a detection of one false CP in every 3 minutes. During the subsequent analysis 
only events that occurred in the narrow windows located at the longitudes of the main pulse 
and the interpulse were selected. These windows were each about lOOyUS wide and occupied 
only 0.6% of the period. Therefore, false detection happened only once per several hours. 

The total number of detected events was about 10000, 2700, 1600, and 1 at the frequen- 
cies of 0.6, 1.4, 2.2, and 8.3 GHz, respectively (see Table 1). At 1.4 and 2.2 GHz sufficiently 
strong CPs were detected simultaneously in several frequency channels, and the number of 
real CPs is 2-3 times less than the total number of detections. The brightest CP was found 
at 1.4 GHz with the peak flux density of 3.7 MJy. Examples of strong giant pulses detected 
simultaneously at 1.4 and 2.2 GHz are shown in Figure 2. 

4. Giant pulse at 8.3 GHz 

Only one giant pulse was detected at the observing frequency of 8.3 GHz at a level of 
30(7 with a peak flux density of about 10 kJy, nearly equal in all four frequency channels. The 
pulse profile is shown in Figure 3. The event was not accompanied by notable emision at any 
of the lower frequency ranges, thus, proving Au/u < 0.5 for the radio spectrum of this pulse. 
On the other hand, at 8.3 GHz the pulse flux density is constant over about 80 MHz. The only 

^ http : //www . jb .man. ac .uk/~ pulsar/ crab .html 

At 600 MHz we did not selected individual GPs (except for the preliminary search). Instead of this, we 
compared the overall statistical properties of signal on emission longitude windows with those on the noise. 
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Fig. 2. Example of strong giant pulses detected simultaneously at 1.4 and 2.2 GHz in several frequency 
channels. The zero value in time corresponds to the sample with the maximum flux density of the giant 
pulse in any of the channels. 

detection did not allow us to perform a statistical analysis of GP shapes, widths, and energies as 
we did for other frequencies (see below). High-frequency observations of giant pulses from the 
Crab pulsar were recently carried out by some of us with the 100-m Effelsberg radio telescope 
at 8.7, 15.1, 22.2 GHz. The first results from the ultra high time resolution data obtained with 
the Tektronix oscilloscope were reported by Jessner et al. (2008). 

5. Giant pulses at 1.4 and 2.2 GHz 

Unlike the high frequency of 8.3 GHz where only one GP was detected, we have found 
1112 and 631 GP events in all four frequency channels at the frequency band of 1.4 and 2.2 GHz, 
respectively (see Table 1). Below we present a statistical analysis of GP energies and widths, 
and a study of radio spectra, shapes, and DM variation. 

5.1. GP width and energy 

Giant pulses vary significantly in their intensity and width, and show a broad variety 
in their shapes as can be seen in Figure 2. Therefore, it is rather difficult to compose an 
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Fig. 3. Giant pulse detected at 8.3 GHz. 

automated routine for GP's shape analysis. In this paper we propose two different techniques 
for such an analysis: the first one is based on a calculation of an integrated flux density 
over varying time boxes in the vicinity of a given GP (this paragraph), and another method 
uses analysis of a signal variance (see paragraph 6.1.1). Let us define the total width of 
a giant pulse (wtot) as a time interval At covering the vicinity of GP which contains 90% 
of the total integrated flux density of a given GP. We will calculate the integrated flux 
density (or pulse energy) as E = 5t(I]"=o('S'i ~ i^)))^ where 6t is sampling time (in our 
case 6t = 0.125/is), Si are individual counts of reconstructed flux density. To exclude the 
noise contribution and calculate the "net" on-pulse energy, one has to subtract the mean 
off-pulse flux density, (5), computed for the off-pulse windowj from flux density value for 
every sample in on-pulse window where the required energy is calculated. In this paper 
we will express pulse energy in units of kJy /is. Following these definitions, we calculated, 
first, the maximum value for required pulse energy separated into two components Eq and 
Eq as E^ = 5t(Er=o('S'max+i - {S))) and E^ = 5t(Er=o('S'max-i - (S))). These quantities are 
calculated in the vicinity of a point of maximum flux density of a GP Smax over n samples in 
positive and negative time directions. We use n = 256, corresponding to maximum expected 
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pulse width At = 2n6t = 64/is. Eq and Eq may differ notably since pulse shape may be 
not symmetrical relative its maximum. With Eq and Eq in hand we calculated trial values 
Ek = ^^(Ei=o('S'max+i - {S))) and E^ = 5t(Ei=o('S'max-i - {S))) for every step k in organized 
loop from k = to n. On every step k we compared the value of Eq with the trial value of 
E'^ , and the value of Eq with the value of E^ . When the corresponding difference \E^ — E^\ 
was less than O.IE^, then the loop was terminated and corresponding time or t~ was 
determined as the edge of a GP. Thus, the total GP width Wtot was defined as Wtot = t'^ —t~. 
Correspondingly, the GP energy Ep was defined as Ep = E^ + Ej^ , and the effective width of 
GP as Wc = Ep/ ^max- Term We is traditionally used in pulsar researches, it describes time inter- 
val where most of pulse energy is concentrated, while Wtot corresponds to the whole pulse extent. 

There is a strong constraint on the usage of this approach to analyse GP's shape. As it 
was already mentioned, the statistics of a detected signal in our case is close to a x^-distribution 
with two degrees of freedom. For such a signal the root-mean-square deviation is equal to the 
mean value, i.e. the modulation index is equal to 1. Thus, the values of E^ will manifest notable 
fluctuations. To keep these fluctuations reasonably small one has to select only sufficiently 
strong pulses. We selected a threshold for the peak flux density for such pulses of 50a for the 
estimated value for maximum pulse width of Wmax = 64 fis. Only 357 and 198 pulses met these 
requirements at 1.4 and 2.2 GHz, respectively. 

Histograms of the total GP width Wtot at 1.4 and 2.2 GHz are shown in Figure 4 (left). 
Both top and bottom histograms are skewed heavily towards shorter total widths. Note that 
only very strong GPs were analysed; a weaker population of GPs may have a broader width 
distribution. Right histogram in Figure 4 represents the dependency between the maximum 
flux density, ^max? and the effective pulse width Wc of a GP. It is clear that the strongest pulses 
tend to have shorter durations. The same tendency was found by Popov & Stappers (2007) 
at a frequency of 1.2 GHz on much broader statistics. We did not analyze the cumulative 
disributions at 1.4 and 2.2 GHz because of rather poor statistics as compared with the data 
used by Popov & Stappers (2007), but we will present such an analysis for our data, obtained 
at 600 MHz in section 6. 

5.2. Instantaneous GP radio spectra 

Given the multifrequency observations with the characteristic frequency setup (see 
Table 1) we were able to analyze instantaneous radio spectra of detected strong GPs both 
over large and short frequency scales. It is important to take into account the effect of in- 
terstellar scintillations when measuring the integrated flux density of a giant pulse in a given 
receiver band. For this we considered results of another experiment that was carried out by 
some of us over the course of our observations. This experiment was the two-station VLBI 
observations between the 64-m Kalyazin radio telescope and 43-m Algonquine radio telescope 
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Fig. 4. Distribution of the total GP width Wtot (left), and dependence between the peak flux density 
'S'max and the effective width of GP Wc (right) at 1.4 and 2.2 GHz. 

(Canada) at 2244 MHz using the S2 recording system in a 16-MHz band (Kondratiev et al. 
2007). In the cross-correlation functions (CCF) of radio spectra of GPs detected at two stations, 
two frequency scales of diffraction modulation of the spectrum at 2244 MHz were found to be 50 
and 450 kHz with equal modulation indices (Kondratiev et al. 2007). Cordes et al. (2004) found 
that in the presence of interstellar scintillations the integrated flux density over the receiver 
band B will be subject to random fluctuations which follow a x^-distribution with OAB/Ai^a 
degrees of freedom, where Aud is the decorrelation bandwidth. In the case of Aud > 0.2B the 
measured integrated flux density is fully modulated by diffractive interstellar scintillations. For 
the VLBI observations having Az/^ = 450 kHz at 2.2 GHz and 5 = 4 MHz, Az/^ > O.IB and 
the modulation of measured integrated flux densities will be still very significant. Following 
the steep power-law frequency dependence of decorrelation bandwidth Az/^, Au^ ~ z/^-^, for the 
Kolmogorov spectrum of electron density irregularities (see, e.g., Lorimer & Kramer, 2004), at 
a frequency of 1.4 GHz, Aud will be about 5 times less than at 2.2 GHz, and the scintillation 
modulations will be considerably smoothed. Therefore, we analysed properties of instantaneous 
GP radio spectra at short frequency scales only for the frequency band of 1.4 GHz. 

The left plot in Figure 5 shows an example of instantaneous GP radio spectra over 20 
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Fig. 5. Instant radio spectra of single GP (left) and their structure functions of time variations for two 
successive observing days (right) at 1.4 GHz. Every curve connecting the same symbols in the left figure 
represent the form of instant radio spectrum of a given GP. The location (shift) of a given spectrum at 
ordinate axis corrcsponts to the arrival time of the GP expressed in minutes from the beginning of the 
observing set. 

minutes of observations at 1.4 GHz. Points presented by the same symbols and connected with 
the same type of a hne represent the pulse energy measured for a given GP and expressed in 
arbitrary units, but with the same scale for an every curve. The every curve is placed at the 
ordinate axis in rough correspondance with arrival time of a given GP. Arrival time is expressed 
in minutes from the beginning of a selected observing interval. The whole figure corresponds 
to one continuum observing set. The variations of presented radio spectra with time seem to 
be random, but one can see that around the 10th minute of observations there are four GP 
spectra similar in appearance. This fact means that interstellar scintillations do still affect the 
measurements. To make a quantitative analysis we calculated a structure function C^(r) of time 
variations of instantaneous GP A structure function is traditionally used in order to classify time 
variations of some quantity Z{t) being investigated. For such a case C^(r) = j^Yl!i=i{Zi — Zi+T-Y . 
Then, C^(r) = with r = 0, and C^(r) = cr^ for r > t with t indicating a characteristic time 
of a variability of Z{t). Here we transform the definition of C^(r) to apply it for the analysis 
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of time variations of 4-point instantaneous spectra F/ as C^(r) = ]g Z^ili Z]/=i (^/ — Y/_^^y. 
Here, M is the number of combinations of radio GP spectra obtained within a given interval 
of time delay r, is the total number of GPs in this test, and / = 1-4 corresponds to one 
of four frequency channels listed in Table 1. The calculated structure functions, C^(r), for 
two successive observing days are shown in the right plot of Figure 5. It is evident that the 
structure functions got saturated at time intervals greater than one minute and decrease rapidly 
at shorter time delays. Thus, we do observe the effect of interstellar scintillations at 1.4 GHz 
that affects our measurements of radio GP spectra. Nevertheless, even at the shortest time 
delay (less than 10 seconds) between instantaneous GP spectra the structure functions still 
have a notable value of about 0.5 that can be due to an intrinsic modulation of GP spectra at 
a short frequency scale of Au/u ^ 0.01. 
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Fig. 6. Instantaneous radio spectra of single GPs in the frequency range of 0.6-2.2 GHz. Every polygonal 
line connecting the same symbols represents the spectrum of a selected giant pulse. 

Examples of instantaneous radio spectra in the frequency range of 0.6-2.2 GHz for several 
GPs are presented in Figure 6. Every polygonal line connecting the same symbols represents 
the spectrum of a selected giant pulse. GP energies at 1.4 and 2.2 GHz were averaged over 
four frequency channels. There were only 27 strong GPs that occurred at all three observing 
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frequencies of 0.6, 1.4, and 2.2 GHz. The mean spectral index ot is equal to —2.6 and —1.8 for 
0.6-1.4 GHz and lA-2.2 GHz, respectively^ Though our estimates are based on poor statistics 
of GPs, the tendency for the spectral index a to flatten at higher frequencies is doubtlessly 
apparent. The same behaviour was reported by Popov et al. (2008) for a frequency range of 
600-1650-4850 MHz. 

5.3. DM variations 

There are 952 GPs out of the total amount of 1112 GP events (or about 85%) found at 
1.4 GHz that were identified with the GPs detected at 600 MHz. Therefore, we used them to 
determine the accurate values of DM by comparing their TOAs at two observing frequencies. 
The results are shown in Figure 7. 
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Fig. 7. Dispersion measure (DM) variations determined from times of arrival (TOA) of GPs at 600 and 
1400 MHz. The left plot shows the measured DM values for different GPs versus the TOA of the giant 
pulse. Open squares show GPs that occurred at the longitude of the main pulse, while filled squares are 
those that occurred at the longitude of the intcrpulsc. The right plot shows the histogram of DM variations. 
The dashed line is the two-Gaussian-componcnt fit with the mean values of 56.7363 and 56.7373 pc cm"'^. 
This difference corresponds to the time shift of about 10 /is at 600 MHz. 



The spectral index a is defined here as E^, oc i>°', where is the integrated fiux density at frequency i>. 
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The left plot of Figure 7 represents the dependence of the measured value of DM of 
GPs versus their time of arrival. Hankins & Eilek (2007) pointed out that at high radio 
frequencies (above 4 GHz) DM values of GPs occurring at the longitudes of the interpulse, 
differ noticeably from DM values of GPs from the longitude of the main pulse. It is clear that 
in our observations at lower frequencies of 0.6 and 1.4 GHz there is no difference between DMs 
for GPs from different longitudes (see Figure 7, left). The right plot of Figure 7 shows the 
histogram of the dispersion measures of GPs. One can notice that the DM distribution is likely 
bimodal and can be well- fitted by two Gaussians with mean values of 56.7363 and 56.7373 pc 
cm~^. The difference ADM = 0.001 ± 0.0001 is estimeted with a good formal accuracy. The 
difference corresponds to the time shift of about 10 fis at the observing frequency of 600 MHz. 
As it follows from the left plot in Figure 7 the splitting can not be attributed to GPs occurred 
at different longitudes (MP and IP). 

5.4- GP profile shapes 

As will be shown below in section 6.3 during our observations the scattering contribution 
was small, and this allowed us to perform the analysis of intrinsic shapes of GPs detected at the 
frequency of 2.2 GHz. Indeed, the expected scattering time Tgcat at this frequency was ^ 200 ns 
as extrapolated from the Tgcat at 600 MHz, which is less than our time resolution of 250 ns 
(doubled sampling interval). 

To study the profiles of individual GPs at the frequency of 2.2 GHz we used the approach 
of calculating the autocorrelation functions (ACFs) of GP profiles, the technique traditionally 
used to study the properties of the microstructure in pulsar radio emission (Hankins 1972; 
Cordes & Hankins 1973; Cordes 1976b). However, to avoid the noise spike at zero lag of 
the ACFs which would mask the presence of any unresolved components in GP profiles, we 
calculated CCFs between the profiles of the same GP recorded in different frequency channels. 
The cross-correlation function -Ri,2(t) between the profiles Ii{t) and l2{t) of the same GP 
detected in two frequency channels is defined as 

N 

RiAr) = [i?i,i(0)i?2,2(0)]-^/'E^iW^2(t + r) , 

t=i 

where t is the time sample within the GP, is the total ON-pulse number of samples, and 
h{t), hit) represent ON-pulse intensities with the OFF-pulse mean level subtracted. In our 
analysis we chose the duration of the ON-pulse window of 125 fis, ot N = 1000 samples. For 
our 4 frequency channels at 2.2 GHz there are C| = 6 different frequency combinations listed 
in the first column of Table 2. 

Due to the effect of interstellar scintillations at this high frequency of 2.2 GHz many 
GPs were not sufficiently simultaneously strong in all frequency channels to provide a notable 
correlation coefficient for every pair of GP profiles. Thus, after the calculation of CCFs for all 
possible combinations of frequency channels for all 631 GP events found at 2.2 GHz, we chose 
only CCFs that have a sufficiently large value of correlation coefficient at zero lag. For an ACF 
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(or CCF) the root-mean-square deviation ((Xacf) near zero lag is equal to (Xacf = 1/VN. In our 
case for = 1000, cTacf ~ 0.03, and for the subsequent analysis we selected only the CCFs with 
the correlation coefficient at zero lag larger than Sa^cf ~ 0.25. There were 121 CCFs selected 
in this way. 

Figure 8 shows an example of selected CCFs calculated for one of the strong CPs. One 
can see the notable correlation between all frequency channels. There are two apparent time 
scales in selected CCFs, namely a short one that occurred as a narrow feature that can be 
attributed to the presence of unresolved spikes in the CP structure, and an extended one with 
a duration of ~ 1 fis. The majority of the selected CCFs had the same appearence. To study 
the extended time scale feature, we smoothed the CCFs over 20 samples and selected only 
those with a correlation coefficient of 0.06 or larger at zero lag. Altogether, 102 of such CCFs 
were selected. The right plot of Figure 8 shows the distribution of both time scales (short 
and extended) measured at the half-level of the maximum of the corresponding broad feature. 
There is a broad jitter of widths of the extended scales in the CCFs from about 1 /is up to 
more than 2 /is with the mean value of about 1.3-1.4 /is. 
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Fig. 8. Examples of cross-correlation functions (CCFs) between profiles in single frequency channels for 
one of the strong CPs at a frequency of 2.2 GHz (left), and the histogram of the half-widths of the two 
time scales detected in the CCFs (right). 
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Table 2. Extra non-dispersive time delay St of the peak in the selected CCFs between GP profiles in different frequency 
channels. These delays represent the residual delays in time of arrival of GPs. 



Frequency combination 


6t 


(MHz) 




2140-2160 


-0.20 ±0.04 


2140-2180 


-0.17±0.04 


2140-2200 


-0.26 ±0.07 


2160-2180 


0.05 ±0.05 


2160-2200 


-0.06 ±0.07 


2180-2200 


-0.01 ±0.08 



In addition to the measurements of the half-widths of the time scales detected in the 
CCFs, we also calculated the time delays of the CCF peaks relative to the expected value for 
the accepted DM. The results are listed in the Table 2. One can see that apparent deviations (at 
a level of > 3a) in GP times of arrival are present only for the frequency channel of 2140 MHz 
(the lowest frequency channel at 2.2 GHz). These deviations measured in a 60-MHz frequency 
range (2140-2200 MHz) at 2.2 GHz are much smaller than the residual dispersion delays of 
about 2 fis found in arrival times of micropulses for the pulsars B0950-I-08 and B1133+16 at 
1.65 GHz over a 16-MHz frequency range (Popov et al. 2002). 

There is also another difference between CCFs computed for GPs in our current analysis 
and ACFs/CCFs calculated for the micropulses (see, e.g., Popov et al. 2002). Parameters of 
the microstructure are explained by the amplitude-modulated noise (AMN) model, developed 
by Rickett (1975). The model explains the pulsar microstructure as an envelope of a complex, 
white Gaussian noise process. Further elaboration of this model was done by Cordes (1976a) 
who postulated the origin of the shot noise as a coherent curvature radiation from bunched 
particles. Cordes (1976b) showed that ACFs/CCFs for micropulses do follow the prediction 
of the AMN model concerning the amplitude of the microstructure component in the average 
ACF/CCF. This amplitude is close to 0.5 at zero lag, when properly corrected for a signal- 
to- noise ratio (S/N) and for a time smoothing. Thus, even shortest micropulses do consist of 
a large number of such shots. However, the CCFs calculated for GPs in our analysis do not 
obey the prediction of the AMN model, and correlation coefficients at zero lag are often more 
than 0.5 level, even without the correction for the finite S/N (see Figure 8). Therefore, shortest 
GP components at 2.2 GHz consist of a small number of individual shots, or even a single 
elementary shot. 
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6. Giant pulses at 600 MHz 



Giant pulses at 600 MHz appear more often than at higher frequencies. In more than 
5 hours of observations we found 997 GPs with peak flux density greater than 30a — 849 in 
main pulse and 148 in interpulse. The highest peak flux density detected was 2 MJy (ITOOcr). 

Owing to unusually low scattering (see section 6.3), we were able to estimate the width 
of the pulses and find the energy distributions separately for pulses of different widths. Below 
we will briefly describe our calculation of energies and widths, and then discuss the results. 

6.1. Data analysis 

It is well known that GPs from Crab pulsar occur in two narrow longitude regions of 
about 6.7° each, called main pulse (MP) and interpulse (IP). Therefore, we analyzed data 
only from three relatively narrow windows: MP, IP and control noise window of the same size 
between MP and IP (with no pulsar emission). 

At 600 MHz the search methods differed from those at other frequencies. Instead of 
selecting the pulses brighter than fixed threshold, we recorded the "width" and "energy" in all 
three windows on every period, regardless of the actual presence of giant pulse. The properties 
of GPs themselves were derived from comparison of the whole sets of "widths" and "energies" 
on MP/IP windows with those of the noise window. We believe such method has two main 
advantages. First of all, it allows us to neglect the unexpected deviations of noise from Gaussian. 
Second, it exempts us from selecting pulses by their peak flux density, which can be randomly 
biased by scattering. Considering all this, the analysis can be extended to weaker pulses. 
6.1.1. Calculation of width and energy 

Giant pulses have complex, jagged form, sometimes with several close components (see 
Figure 9). Here we considered another technique of CP's width measurements compared with 
that proposed in paragraph 5.1. To find the width of CP, we used the fact that every giant 
pulse, regardless of its shape and peak flux density, is, in fact, a region with the increased 
variance of the signal S. Thus, on every MP/IP/noise window we computed a sequence of 
variances, af over small floating window: 

N 
k=l 

The size of the floating window, A^, was comparable to the typical distance between the 
components (namely 200 samples, i.e. 25fis). When such window covered CP, the variance 
rapidly increased. Since the size of the window was about the typical distance between the 
components of GP, the raise and decay were smooth and it was easy to measure the width of 
such outburst which we designated as W^2. It was measured on the level of 5% of the peak 

i 

maximum (see Figure 9). To obtain the width of GP (Wgp) we subtracted the size of the 
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window from W„2: 

Wgp = W,2-N . 

i 




N Wqp 



Running variance, a-^ 



-.- . -r I 

200 400 600 800 1000 1200 1400 

Number of counts, i 

Fig. 9. Example of single pulse detected at 600 MHz, illustrating the technique of pulse width definition. 
Solid line represents the pulse intensity, and dash line shows the runninig variance. See text for more 
detailed explanation. 

After that, the signal on corresponding MP/IP or noise window was smoothed with the 
value of the width. The maximum on the smoothed array was regarded as the total energy Etot 
of an event, and obtained Wgp is understood as We- Then we converted derived We, expressed 
in samples, to microseconds with 6t = 0.125/is, and we scaled Etot to kJy ■ fis. 

As was mentioned before, such procedure was performed in every period in three different 
longitude windows. On noise window W^2 was usually less than the size of floating window 

i 

N, resulting in negative "GP" widths. To unify the smoothing procedure we assigned all such 
negative widths the value of 1. 
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Fig. 10. Histogram of total widths (left), and cumulative energy distribution for GPs with widths of 16-24 
fis at 600 MHz (right). Dash line in the histogram represents GPs with peak flux densities Speak > 50(7 
(the strongest GPs), while solid line is for weaker GPs. 

6.2. GP widths and energies 

As can be seen on Figure 10, the strongest pulses are, on average, narrow, and their 
width is comparable to the scattering width. The dimmer pulses are usually broader. The 
same tendency was observed at 1.2 GHz (see section 5). 

The energy distributions for GPs both from the main pulse and the interpulse are well 
fitted by a power-law, with the indices depending on the range of widths selected (Table 3). 
However, for MP GPs wider than 8 /is, the energy distribution flattens at low energies (see 
Figure 10), while energy distribution for IP GPs of the same width and energy remains the 
same. Thus, the cumulative distribution for GPs from main pulse wider than 8 /is was fitted 
with 2 power-law functions: 



N{E > Eo 



Eq^ < -E'break! 

Eq'^ , Eq > -Ebrcak 



Such flattening was previously detected at other frequencies (Bhat et al. 2008; Popov & 
Stappers 2007 and references therein). The break energy for overall distribution, 35 kJy-/is, is 
well consistent with the empiric formula for the break energy in dependance on the observed 
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Table 3. The power-law indices a of the cumulative energy distributions of GPs at the frequency of 600 MHz for different 
GP widths. Energies of GPs from the longitude of the main pulse have different power-law distributions for low-energy and 
high-energy GPs with the clear break energy -Ebrcak (see text) between them. We list both indices ai and 02 for this case. 



iotai (jr WlQin, 


Intcrpulsc 




Main Pulse 




a 


Oil 


a2 


-^break i 

Jy-pis 


<8 


-1.5 




1.4 




8-16 


-1.9 


-1.2 


-1.9 


40 


16-24 


-2.4 


-1.3 


-2.4 


45 


> 24 


-3.2 


-1.0 


-3.8 


55 


all 


-2.0 


-0.9 


-2.2 


35 



frequency, derived by Popov & Stappers (2007): 

-E'break = ^'^ , 

where E is in kJy -/is and u in GHz. 
6.3. GP scattering 

At frequencies below ~ 1 GHz pulse shapes are dominated by the scattering broadening 
which decreases with increasing frequency for the Crab pulsar as Tgcat oc /~^'^ (Popov et al. 
2006a). However, the amount of scattering varies significantly with time. For example, Popov 
et al. (2006b) found Tscat to be equal to 45 /is at the frequency of 600 MHz in their observations 
with the Kalyazin radio telescope in November, 2003. The measured value of Tscat was estimated 
by fitting the exponential function exp(— t/rscat) to the tail of the average GP profile. In our 
current observations in July, 2005 the scattering broadening was 5-10 times smaller, and to 
approximate the average GP profile tail we used two exponents with the Ts^at of 4.5 and 17.5 /is, 
with short-scale component being 4 times larger in amplitude. The mean GP profile was folded 
by summing 1436 GPs stronger than 50 kJy in peak flux density. They were preliminary 
averaged by 10 samples (1.25 /xs), and aligned by their abrupt leading edges. The tail portion 
of the average GP profile at 0.6 GHz and its fit are presented in Figure 11. Smaller scattering 
contribution favoured our analysis of intrinsic shapes of GPs detected at 2.2 GHz (Section 5.4). 

7. Conclusions 

Accomplished study of giant pulses from the Crab pulsar from our simultaneous mul- 
tifrequency observations with the Kalyazin radio telescope allowed us to draw a number of 
conclusions about the properties of Crab GPs: 

1. Single GPs were often observed simultaneously in the three frequency bands of 0.6, 1.4, 
and 2.2 GHz, confirming the broad band nature of the GP emission. The spectral index a 
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Fig. 11. The tail portion of the average GP profile at 600 MHz folded by summing 1436 GPs 
stronger than 50 kJy in peak flux density. The dashed line represents the fit by two exponents 
S = 8.1exp(-t/Tscat,i) + l-8exp(-t/Tscat,2) with the scattering times r^cat,! = 4.5 fis and rscat,2 = 17.5 ^s. 

has a tendency to flatten towards higher frequencies (ao.6-1.4 = —2.6 and ai.4_2.2 = —1-8) 
as was also reported before by Popov et al. (2008). 

2. Instant radio spectra of individual GPs do not follow the single power law in the frequency 
bands of 0.6, 1.4, and 2.2 GHz. Many GPs observed both at 0.6 and 2.2 GHz were not 
detected at 1.4 GHz, thus proving a notable spectra modulation at a frequency scale of 
Au/u ^ 0.5. At the frequency of 1.4 GHz in addition to diffraction spectra distortions 
caused by the wave propagation through inhomogeneties of the interstellar plasma, we 
also distinguished the small-scale modulations in radio spectra at a level of Au/u ^ 0.01 
intrinsic to the pulsar. 

3. GP profiles at the frequency of 2.2 GHz can be presented by unresolved spikes grouped 
together at a time scale of about 1 /is. These unresolved components are well-correlated 
over 60-MHz bandwidth (Az//z/ ^ 0.03). For comparison, Hankins (2000) did not find 
the correlation for the unresolved structure of GPs between 4535 and 4985 MHz (i.e., 
Au/u^O.l). 

4. Cross-correlation functions of GP profiles in different frequency channels do not obey 
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the prediction of the AMN model (Rickett 1975) for the microstructure of the pulsar 
radio emission. Thus, unresolved components represent a small number or even the single 
elementary emitter. 

5. The extra non- dispersive time delays of about 0.2 fis found in arrival times of unresolved 
components of GP structure at 2.2 GHz are not notable and 10 times smaller than the 
residual dispersion delays of about 2 /is found in TOAs of micropulses for the pulsars 
B0950+08 and B1133+16 at 1.65 GHz (Popov et al. 2002). 

6. Strongest GPs tend to have shorter durations. Same results were reported by Popov & 
Stappers (2007) and by Bhat et al. (2008) at frequencies of 1.2-1.4 GHz. 

7. The distribution of dispersion measures of GPs determined from their TOAs at 0.6 and 
1.4 GHz is likely bimodal and can be well-fitted by two Gaussians with mean values of 
56.7363 and 56.7373 pc cm~^. 

8. The cumulative distribution of GP energies at 600 MHz depends on the effective width 
of GPs. This distribution obtained only for GPs at the longitude of the main pulse, has 
the break and can be represented by two power-law functions with indices of -2.0 - -0.9 
at high and low energies respectively. The value of the energy -Ebreak, where this break 
occurred, is equal to 35 kJy-/is, and it fits to the empiric formula inbreak (kJy ■ /is) ~ 7i'~^-'^ 
with u in GHz derived by Popov & Stappers (2007). 
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